The infrared behaviour of the static potential in perturbative QCD 
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I. INTRODUCTION 

The static singlet potential between a heavy quark 
and an antiquark is an object of considerable interest 
both theoretically and phenomenologically. Strictly 
speaking, being the potential a dynamical quantity, 
it can be defined only in a context where the full 
dynamics of the system is consistently taken into 
account. For instance in the proper Schrodinger 
equation. Nevertheless, the static singlet potential is 
usually defined as the logarithm of the static Wilson 
loop divided by the interaction time in the limit of 
infinite interaction time. Moreover, the above defini- 
tion, in QCD, suffers from infrared (IR) divergences 
when computed at finite orders of perturbation the- 
ory HI . This has to do with the non- Abelian nature 
of QCD, which allows massless gluons to self-interact 
at arbitrarily small energy scales. 

The problem was posed in the past jjj of whether 
the static singlet potential could be defined in some 
way at any order in perturbation theory. The lead- 
ing IR singularities in the static Wilson loop were 
found at a relative order a 3 in Q. These singulari- 
ties can indeed be regulated upon resummation of a 
certain class of diagrams which give rise to a dynam- 
ical cut-off provided by the difference between the 
singlet and the octet potentials. However, the ex- 
traction of the static singlet quark-antiquark (Q-Q) 
potential from this resummation has been regarded 
as suspect by other authors [Q, who noticed that 
the difference between the singlet and the octet po- 
tential provides a dynamical scale which, at least in 
perturbation theory, is of the same order of the ki- 
netic energy for quarks of large but finite mass. A 
similar situation occurs in QED for the Lamb shift, 
where one of the contributions has to be interpreted 
as a non-potential correction to the hydrogen atom 
energy levels. The novel feature in QCD is that such 
an effect appears to be related to the definition of the 



static potential, while in QED it is of order 1/m 2 . 
No quantitative procedure has been developed so far 
to deal with this problem in QCD. 

The recent complete calculation of the static po- 
tential at two loops H and its foreseen applications 
to top pair production at the Next Linear Collider 
|Q brings this problem to the edge of experimentally 
testable physics. We feel that a clear resolution is 
urgently needed. 

It is the aim of this paper to fully clarify this 
problem by providing a quantitative framework 
where the different scales that characterize a non- 
relativistic system can be properly taken into ac- 
count. This is achieved by using an effective field 
theory approach. 

We will assume that mv 3> Aqcd, wi and v be- 
ing the heavy-quark mass and velocity respectively 
(v <C 1). Then we will systematically integrate out 
the hard (~ m) and soft (~ mv) degrees of freedom 
by performing in both cases a perturbative matching 
to suitable effective field theories. The integration 
of the hard scale gives rise to non-relativistic QCD 
(NRQCD) ||, whereas the integration of the soft 
scale gives rise to potential NRQCD (pNRQCD) §]. 
In the latter effective theory only ultrasoft (US) de- 
grees of freedom are left, which in this letter means 
the ones with energy much smaller than mv. The 
static potential can be understood as a matching 
coefficient in pNRQCD. 

The main result of this work is to state rigorously 
what in perturbative QCD has to be understood as 
the Q-Q static potential, namely the relevant ob- 
ject for the dynamics of Q-Q pairs with large but 
finite mass. This does not simply coincide with the 
static Wilson loop as usually computed, since this 
turns out to contain US contributions as well. As 
a consequence the static potential manifests, at the 
relative order a 3 , an explicit dependence on the cut- 
off of the effective theory and has to be understood 
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as a matching coefficient in pNRQCD. The leading 
cut-off dependence is evaluated explicitly. 



II. THEORETICAL FRAMEWORK 

After integrating out the hard scale (~ m) from 
QCD, we are left with NRQCD. For our purposes 
it is sufficient to work at the lowest order in the 
NRQCD Lagrangian, namely, 



'4 » v 



where ip is the Pauli spinor field that annihilates the 
fermion and \ is the Pauli spinor field that creates 
the antifcrmion; iD° = ido — gA° and iT) = iV + <?A. 

Integrating out also the soft scale, mv, from (Q) we 
are left with an effective theory (pNRQCD) where 
only US degrees of freedom remain dynamical. The 
surviving fields are the Q-Q states (with US en- 
ergy) and the US gluons. The Q-Q states can be 
decomposed into a singlet (S) and an octet (O) un- 
der colour transformation. The relative coordinate 
r = xi — X2, whose typical size is the inverse of 
the soft scale, is explicit and can be considered as 
small with respect to the remaining (US) dynamical 
lengths in the system. Hence the gluon fields can be 
systematically expanded in r (multipole expansion) . 
Therefore the pNRQCD Lagrangian is constructed 
not only order by order in 1/m, but also order by 
order in r. As a typical feature of an effective the- 
ory, all the non-analytic behaviour in r is encoded in 
the matching coefficients, which can be interpreted 
as potential-like terms. 

The most general pNRQCD Lagrangian density 
that can be constructed with these fields and that is 
compatible with the symmetries of NRQCD is given 
at the leading order in the multipole expansion by: 
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pNRQCD 



= Tr\sUid -?--V s (r) + 
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+0t [iDo-^-V (r) + ...)Oj 
+gFA(r)Tr{O t r-ES + St r -EO} 
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TrjO^-EO + OtOr-E}, (2) 



where R = (x a + x 2 )/2, S = S(r,R,t) and O = 
0(r,R, t) are the singlet and octet wave functions 
respectively. All the gauge fields in Eq. (Q) are eval- 
uated in R and t. In particular E = E(R, t) and 



iD O = id O - g[Ao{R,t),0]. V s and V a are the 
singlet and octet heavy Q-Q static potential respec- 
tively. Higher-order potentials in the 1 jm expansion 
and the centre-of-mass kinetic terms are irrelevant 
here and are neglected. We define 



V s (r) = -C } 



a Vs (r) 



(3) 



Va and Vb are the matching coefficients associ- 
ated in the Lagrangian (|^) to the leading correc- 
tions in the multipole expansion. Both the poten- 
tials and the coefficients Va and Vb have to be 
determined by matching pNRQCD with NRQCD 
at a scale fj, smaller than mv and larger than the 
US scales. Since, in particular, [i is larger than 
Aqcd the matching can be done perturbatively. At 
the lowest order in the coupling constant we get 
ay 3 — ctv — a s , Va = Vb = 1- In order to have the 
proper free-field normalization in the colour space 
we define 

1 T a 
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where Tp = 1/2. 



III. MATCHING 

In this section we discuss how to perform the 
matching between NRQCD and pNRQCD. In par- 
ticular we will concentrate on the effects produced 
by the leading corrections coming from the multipole 
expansion. 

The matching is in general done by comparing 2- 
fermion Green functions (plus external gluons at the 
ultrasoft scale) in NRQCD and pNRQCD, order by 
order in 1/m and order by order in the multipole 
expansion. If the soft scale is in the perturbative 
region of QCD (i.e. larger than Aqcd), this can be 
done explicitly order by order in the coupling con- 
stant. If not, one can still perform the matching by 
subtracting the calculation in pNRQCD to the de- 
sired order of accuracy. The remaining term only 
contains the soft scale (up to US higher-order cor- 
rections) and goes in the pNRQCD Lagrangian as a 
new potential term. 

The matching can be done once the interpolating 
fields for S and O a have been identified in NRQCD. 
The former need to have the same quantum numbers 
and the same transformation properties as the latter. 
We choose the following definitions for the singlet 
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(x 2 , i)#x 2 , x X) t)V>(xi, t) = z¥ 2 (r)S(R, r, t) 
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and for the octet 



X f (x 2 , i)0(x 2 , R, <)T>(R, x 1( t)V(x lf t) 
= Z 1 J 2 (r)O a (R,r,t), 

where 

„ y , X , i ).P« P { !! ,/'„ (y - x ,.A, x - s , x - y ), ! ^ 

From the normalization condition (JI]) it follows, at 
the tree level, that Z s = N c and Z Q = Tp. 

In order to get the singlet potential, we choose the 
following Green function: 

1= {x 2 )cj){x2,xi)^{x l ) 

xtf (yi)ct>(yi, |ft)x(lte)|0). (5) 

In NRQCD we obtain, at order (l/m)°: 

7 = <5 3 ( Xl -yi)(5 3 (x 2 -y 2 )(^ D ), (6) 

where Wn is the rectangular Wilson loop with 
edges x x = (T/2,r/2), z 2 = (T/2,-r/2), yi = 
(-T/2,r/2) and y 2 = (-T/2,-r/2). The symbol 
( ) means the average over the gauge fields. 



W(r x T) 



NRQCD pNRQCD 

FIG. 1. The matching of the static potential. On the 
right side are the pNRQCD fields: simple lines are sin- 
glet propagator, double lines are octet propagators, cir- 
cled-crosses are the singlet-octet vertices of Eq. (^) and 
the wavy line is the US gluon propagator. 

In pNRQCD we obtain at order (l/m)° and at the 
next-to-leading order in the multipole expansion 

I = Z s (r)S 3 ( Xl - yi )£ 3 (x 2 - y 2 )e- lT ^« 
x ( 1 _ 2-Vl(r) / dt / dt'e-^ ^ v °- y ^ 

N r . J-T/2 J-T/2 



:{Tr{v-E(t)4,(t,t')vE(t')^t',t)}) 



(7) 



Fields with only temporal argument are evaluated 
in the centre-of-mass coordinate. 

Comparing Eqs. (||) and (Q), one gets at the next- 
to-leading order in the multipole expansion the sin- 
glet wave-function normalization Z s and the singlet 



static potential V s (see Fig. Q). Va{t) and V (r) 
must have been previously obtained from the match- 
ing of suitable operators. Since here we shall only 
need the tree-level results, we postpone the general 
discussion to We shall further concentrate on 
the singlet potential while the wave function nor- 
malization will also be discussed in |Q . The explicit 
calculation can be done in several (equivalent) ways: 

i) By matching order by order in a s 0. Both 
expressions (||) and @) are IR-divergent (Eq. (Q) 
is also UV- divergent). Regulating both expressions 
in dimensional regularization the calculation in pN- 
RQCD gives zero (there is no scale), while the cal- 
culation of the Wilson loop shows up an explicit de- 
pendence on the infrared regulator /i via a typical 
In fir term. 

ii) By keeping, without expanding, in (0) the ex- 
ponentials in V and V s . The scale V — V s , which 
appears in this way, regulates the IR divergences 
in Eq. ([?]). Because of this scale, the calcula- 
tion in pNRQCD now gives a non-zero contribu- 
tion proportional to the UV cut-off /i of the the- 
ory (~ ln(U — V s )/fi). This calculation would be 
sufficient to extract the leading IR divergences of 
the static potential, in the same way as it is suf- 
ficient to evaluate the leading logarithm in the ef- 
fective field theory in order to know the running of 
a matching coefficient. However, in order to carry 
out the matching consistently, we have to evaluate 
(j^) in a way that exactly corresponds to keeping the 
exponentials in (Q). For the IR-singular terms in 
NRQCD, this way is nothing but the resummation of 
the diagrams depicted in Fig. ^, as it was indicated 
in ffl. This procedure is automatically adopted in 
any attempt to extract the (non-perturbative) static 
potential from a non-perturbative evaluation of the 
Wilson loop (e.g. in lattice calculations). 
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FIG. 2. Graphs contributing to 5{Wa). Dashed lines 
represent Coulomb exchanges. 

In order to make the comparison with Q easier, 
we will follow here this second approach. We will 
show the cancellation of the ln(U — V s ) terms ex- 
plicitly. Evaluating expressions (^|) and (0) with the 
prescription ii) we get, for the static potential at the 
next-to-leading order in the multipole expansion: 
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lim -\n(W a ) 



+c F — lr <y -v a ) 3 In 
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(8) 



The dependence on ln(V D — U s ) in the second line is 
canceled in the first line by the contribution to the 
Wilson loop (5(Wu)) coming from the graphs shown 
in Fig. H whose qualitative features were considered 
in . Indeed, the exact calculation done in this 
work gives 



lim LlnS(Wn) = -C F Cl^-(V 



V, 



xln((U -U s )V) . (9) 

Finally, the result given in Eq. (||) allows us to 
write the complete expression for the static poten- 
tial, Eq. ([}]), up to the relative order o^ln/ir in 
coordinate space (a s is in the MS scheme): 



a Vs (r, fi) = a s (r) < 1 + (ai + 2-f E (3 ) 



+ 



7b (4ai/3 + 2ft) + ( — + 4 7 ^ 
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(10) 



where [3 n arc the coefficients of the beta function, a\ 
was calculated in Ref. ||] and a 2 in Ref. ||] (see Q for 
notation). We emphasize that this new contribution 
to the static potential would be zero in QED. 



separation between soft and US contributions is not 
an artificial trick but a necessary procedure if one 
wants to use the static potential in a Schrodinger- 
like equation in order to study the dynamics of Q-Q 
states of large but finite mass. In that equation the 
kinetic term of the Q-Q system is US and so is the 
energy. Since the US gluons interact with the Q-Q 
system, their dynamics is sensitive to the energies of 
the (non-static) system and hence it is not correct 
to include them in the static potential. When cal- 
culating a physical observable the fi dependence in 
(10) must cancel against /i-dependent contributions 
coming from the US gluons. 

Finally it is worth mentioning that the static po- 
tential suffers from IR renormalons ambiguities with 
the following structure 

5V S ~C + C 2 r 2 + ... 

The constant C is known to be cancelled by the IR 
pole mass renormalon [jl0| while the IR C2 renor- 
malon gets cancelled by the IR renormalon existing 
in the second term of Eq. (||) H . 
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IV. CONCLUSIONS 

Eq. (]Io| ) is the main result of this letter. It states 
that ay , defined through the static potential (see 
Eq. (|j)), is not a short distance quantity as a s 
(in the MS scheme), since it depends on the IR be- 
haviour of the theory. It can be better understood 
as a matching coefficient (an analogous situation oc- 
curs for the pole mass). Moreover Eq. ( [To| ) gives the 
analytic value of the coefficient of the leading In \ir 
correction, which arises at relative order a^, i.e. im- 
mediately after the known two-loop corrections. The 
situation for the static octet potential seems to be 
similar [Q. 

The evaluated terms clarify the long-standing is- 
sue of how the perturbative static potential should 
be defined at higher order in the perturbative se- 
ries. Eq. (||) explicitly shows that the static poten- 
tial does not coincide with the static Wilson loop as 
usually computed. It should be emphasized that the 
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